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We treat the ground state, elementary excitations, and neutron scattering cross section for a 
system of trimers consisting of three tightly bound spins 1/2 on a distorted Kagome lattice, subject 
to isotropic nearest neighbor (usually antiferromagnetic) Heisenberg interactions. The interactions 
between trimers are assumed to be weak compared to the intra trimer interactions. We compare 
the spin-wave excitation spectrum of trimers with that obtained from standard spin-wave theory 
and attribute the differences at low energy to the fact that the trimer formulation includes exactly 
the effects of intra-trimer zero point motion. Application to existing systems is briefly discussed. 
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I. INTRODUCTION 

Frustrated antiferromagnetic systems have re- 
ceived enormous attention in recent years J^— One 
limit which has attracted less attention is that 
when the frustration is removed by the forma- 
tion of strongly coupled three-spin units called spin 
trimers 4^— Early experiments and calculations were 
performed for high (S = 5/2) spin states of Fe 3+ 
and Mn 2+ by Falk et al.jS and Furrer and Giidel. 5 
For S = 1/2 systems much work has been focused 
on chain-like systems consisting of trimers of Cu 
ions4&i£ Other configurations of trimers were stud- 
ied by Qiu et al.j£ and Podlesnyak et al£ In these 
works the interactions between trimers were very 
weak, so that the energy of the localized excitations 
appeared not to depend on wave vector. In that 
case, information on the nature of the excited states 
of the trimers was obtained by monitoring the de- 
pendence of the magnitude of the inelastic scatter- 
ing cross section on wave vector. In contrast, here we 
will consider a system of interacting spin 1/2 trimers 
where the excitations have a significant dependence 
on wave vector. We implement perturbation theory 
by introducing operators which create or destroy the 
the exact excited states of isolated trimers. In the 
limit when the inter-trimer interactions vanish, our 
calculation reduces to those of Refs. Hand0. 

The system of trimers of spins 1/2 we consider is 
specified by Fig. [T] where we show the covering of 
a distorted Kagome lattice by trimers. The lattice 
has the connectivity of a Kagome lattice, but lacks 
its three-fold symmetry, so that the nearest neighbor 
isotropic exchange interactions assume three values 
J, j, and k, of which J is assumed to be dominant. 
This model may be an appropriate one for the dis- 
torted Kagome system Cu2 (OD^CL 11 ' 12 Even if this 
system is not an ideal representative of the model 



we introduce below, our results may stimulate the 
search for better realizations of our model. The aim 
of this paper is to develop a calculation which is 
correct to leading order in j / J and kj J and to com- 
pare results obtained in this approximation to stan- 
dard spin- wave theory, based on the Neel state which 
treats all the exchange interactions on an equal foot- 
ing. We find that there is a one-to-one mapping con- 
necting the lowest energy manifold of excitations in 
the two approaches and that the differences in ener- 
gies can be understood in terms of the differing way 
quantum zero-point motion is treated in the two ap- 
proaches. At higher energy the comparison is more 
complicated. In the trimer approach one does have 
the higher energy transverse spin waves of the Neel 
state. But in addition, some of the higher energy 
trimer excitations correspond to bound states of two 
or more Neel-state spin excitations. The trimer ap- 
proach is clearly superior when the intertrimer inter- 
actions are perturbative, as we assume in this paper. 

Briefly, this paper is organized as follows. In Sec. 
II we give a qualitative overview of the calculation 
in which the intertrimer interactions j and k are 
treated perturbatively with respect to the strong in- 
tratrimer interaction J. In Sec. Ill we show that the 
low energy manifold of spin waves can be mapped 
onto the usual manifold of spin waves, but with an 
effective trimer-trimer interaction playing the role of 
the usual spin-spin interaction. Here and in succeed- 
ing sections we treat the two cases when a) the net 
spins of adjacent trimers are coupled antiferromag- 
netically and b) the net spins of adjacent trimers are 
coupled ferromagnetically. In Sec. IV we consider 
the exciton spectrum in which trimers are promoted 
into their nearly localized excited states. In Sec. 
V we present results of standard spin-wave calcula- 
tions based on the Neel state in which all spins in the 
ground state have S z = 1/2 or S z = —1/2. In Sec 
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FIG. 1: (Color online) A distorted Kagome lattice with 
three isotropic nearest neighbor Heisenberg interactions; 
J (red), j (blue) and k (green). We assume J is antifer- 
romagnetic and much larger than j and k, yielding spin 
trimers (some of which are shown as ellipses) which con- 
sist of three spins connected by two large interactions, 
J. The dashed rectangle is the unit cell containing two 
trimers, A and B. The upper right inset shows the four 
nearest neighbor vectors of the trimer lattice as given 
in Eq. (|14[) . Here s is the nearest neighbor separation 
between spins on the Kagome lattice. The labeling of 
the three sites within a trimer is a, b, c in the order of 
decreasing z coordinate, as shown for a trimer in the left 
bottom corner of the unit cell. 



VI we compare the results the spin-wave and per- 
turbative approaches give for the elastic diffraction 
pattern. We attribute the differences in results to 
the differences in how quantum zero-point motion is 
treated in the two approaches. In Sec. VII we con- 
sider the inelastic neutron scattering cross section 
from the entire spectrum of trimer excitations. Our 
results are summarized and briefly discussed in Sec. 
VIII. 



II. OVERVIEW 

In the magnetically disordered phase of 
Cu(OD)3Cl (which we take as the exemplar of 
our trimer model) the unit cell shown in Fig. [JJ 
contains six Cu spin sites. In Fig. [2] we show the 
phase diagram of the trimer model as a function 
of the temperature T when J is much larger than 
either j or k. When T is large compared to J the 
spins are essentially uncorrelated. As T is reduced 
to become comparable to J, one passes through a 
regime in which the correlations within spin trimers 
become well developed. In Fig. [2 this regime is 
labeled "trimer melting." Below this regime the 
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FIG. 2: (Color online) The phase diagram of the trimer 
system as a function of temperature T, as discussed 
in the text. Long range magnetic order occurs at T c . 
Trimer formation occurs over the regime for which T is 
of order J. 



average spin of the middle site of the trimer is 
oppositely oriented to those of the end sites of the 
trimer. [nlGJ] However, as long as T > T c , the 
spin correlation function between different trimers 
decays rapidly as a function of their separation. 
When T is reduced so as to be comparable to j 
and/or k, one passes through a phase transition 
(at T = T c ) below which one has long range spin 
order. As we discuss below, depending on how j 
and k compare, the adjacent trimers can either be 
organized ferromagnetically or antiferromagncti- 
cally. In either case, the magnetic ordering occurs 
at zero wave vector. In other words, the magnetic 
and paramagnetic unit cells are identical, each 
containing two trimer units. As we shall see, when 
T -C T c the elementary excitations are identical to 
spin waves in the usual magnetic systems. 

In contrast and as will become apparent, the 
higher energy trimer excitons are qualitatively dif- 
ferent from the higher energy spin wave relative to 
the Neel ground state. To obtain a close correspon- 
dence between the two approaches one should con- 
sider trimers consisting of three large S spins. In 
that case, one should pass continuously between the 
trimer and Neel limits as the ratio of j or k to J is 
varied. 

The Hamiltonian for the system of spins 1/2 which 
we treat is written as 



H 



Si ■ 



(i) 



where (ij) indicates that the sum is over pairs of 
nearest neighbors on the Kagome lattice. Here we 
neglect exchange anisotropy, in particular we do 
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not include the Dzialoshinskii-Moriya [13.14J inter- 
action, which can be the dominant anisotropic in- 
teraction between spins. The values of the J's 
are defined in Fig. Q] where the intra-trimer inter- 
action J is assumed to be dominant. We will work 
to leading order in j/J or fc/J which are assumed 
to be of order x « 1. Thus the expansion parame- 
ter x characterizes the ratio of inter-trimer to intra- 
trimer interactions. When inter-trimer interactions 
are turned on, the spectrum of discrete energy lev- 
els of isolated trimers gets broadened into a band 
of wavelike excitations, just as happens for atomic 
energy levels when placed in a solid. 

For this calculation we obviously need the exact 
eigenfunctions and eigenvalues of the trimer Hamil- 
tonian 

H T = JS a ■ S b + JS b ■ S c , (2) 

where the spins within a trimer are labeled as in 
Fig. [TJ The total spin S is a good quantum number 
and assumes the values 3/2 and 1/2. The four states 
5 = 3/2 are degenerate eigenstates of Ht with eigen- 
value J/2. The remaining four eigenstates form two 
5 = 1/2 doublets. The eigenstates and eigenvalues 
of Ht are listed in Table I. In the next section we 
consider the ground state manifold and in the fol- 
lowing sections we consider excitations to the higher 
manifolds centered at energy J and 3 J/2 above the 
ground state. 

Before starting the calculation we should dis- 
cuss when the trimer limit we consider is appro- 
priate. First of all, our results show the obvious 
fact that when the trimers interact with one an- 
other, the single-trimer energy levels get broadened 
into a band. Clearly, a condition for treating iso- 
lated trimers as a starting point, would be that this 
broadening is small enough that the bands are sep- 
arated and qualitatively retain their identity from 
the noninteracting limit. But additionally, in view 
of the fact that the trimers will be shown to act as 
spin 1/2's, one could question whether this calcula- 
tion improves the treatment of quantum zero-point 
which can be severe for S — 1/2. The following qual- 
itative estimate indicates why the trimer calculation 
can be useful. Let us consider excitation relative to 
the Neel state in which spins are aligned along the z- 
axis. The perturbation which creates zero-point mo- 
tion comes from terms like JijS-(i)S+(J)/2, where 
the subscript labels the Cartesian component of spin 
and the largest such terms are those for which sites i 
and j are inside the same trimer. This perturbation 
V connects the ground state to a state with excita- 
tion energy E = 2zJS, where z, the number of near- 
est neighbors should be taken to be 1 or 2 because 
for each site there are only 1 or 2 strongly coupled 
neighbors. Thus V/E 1/2. In contrast, when this 



type of calculation is repeated for the trimer state z 
is now 4, the number of trimer-trimer nearest neigh- 
bors. Also, perturbative corrections to a system of 
isolated trimers are of order V/E = j/J, where j is 
one of the inter-trimer interactions. So zero point 
corrections are less important for the trimer analog 
of the Neel state than for the usual Neel state in the 
limit when j/J is small. 

III. GROUND-STATE EXCITATIONS 

We first consider the 2 Ar -fold degenerate mani- 
fold of ./V trimers when intertrimer interactions are 
turned off, so that each trimer has energy — J. To 
implement degenerate perturbation theory when in- 
tertrimer interactions are turned on, it is convenient 
to map this manifold of states onto the 2^ states 
associated with a system of N pseudospin 1/2 op- 
erators, such that the pseudospin operator of each 
trimer is simply the total ground state spin opera- 
tors S of that trimer. For the trimer at position R 
we denote this pseudospin operator as <r(R). Then, 
any operator within the ground manifold can be ex- 
pressed in terms of products of one or more <r(R). 
We then use the wavefunctions in Table I to express 
matrix elements of spin operators for individual sites 
within the trimer at R to er(R). For this purpose 
we label the three spins within a trimer as a, b, and 
c as in Fig. [T] Using Table U we note that for an A 
trimer (for which cr z = 1/2), the expectation value 
of the z-component of the fcth spin in the ground 
state of the trimer denoted S z (k) (where k = a,b, c) 
is 

S z (a) = 1/3, £,(&) = -1/6, S z (c) = 1/3 .(3) 

This result reflects the fact that the central spin par- 
takes of spin fluctuations with its two neighbors in- 
side the trimer whereas an end spin of the trimer 
has only one neighbor with which to fluctuate. Be- 
low we will discuss the experimental consequences of 
this result. In fact, the Wigner-Eckart theorem^! in- 
dicates that we have, as an operator equality within 
the ground manifold, that 

S(a;R) = 2er(R)/3 , 

S(6;R) = -<r(R)/3 , 

S(c;R) = 2er(R)/3 , (4) 

where S(fc; R) is the operator for the fcth spin in the 
trimer whose center is at R and, as we have said, the 
pseudo-spin operator is identified as the total spin of 
the trimer: 

c 

<r(R) = ^S(fc;R) . (5) 

k—a 
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TABLE I: Eigenvectors tp n and eigenvalues X n of Ht- 
The states specified by three vertical arrows give the 
values of S z for spins a, b, and c (reading from left to 
right), as shown in Fig. [1] The index n is only used to 
label excited states. 
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These equalities make it a trivial matter to write the 
inter-trimer interactions in terms of the er's. So we 
see, even without calculation, that the low-energy 
spectrum of the trimer system is identical to that 
of a system in which each trimer is replaced by an 
ordinary spin 1/2. 

We now consider the ground state and elementary 
excitations of the system when weak interactions be- 
tween trimers are included. We will assume that 
all end-to-end exchange interactions between near- 
est neighbor trimers assume a common value k and 
those between the end of one trimer and the cen- 
ter of its nearest neighbor assume a common value 
j as shown in Fig. 1. This symmetry we have im- 
posed makes the calculations algebraiclly simple. If 
the intra-trimer and inter-trimer interactions have 
no special symmetry, the calculations becomes alge- 
braically more complicated but are conceptually no 
more difficult. So here we give results only for the 
model of Fig. [TJ 

We now construct the effective Hamiltonian within 
the ground state manifold. Consider the interac- 
tion V(A, B) between trimers A and B. We use the 
Wigner-Eckart theorem to express the spin operators 
in terms of the pseudo or total spin of the trimer, as 
done in Eq. Then one sees that V(A, B) within 
the ground manifold is given by 

V(A,B) = [<r(A)-tr(B)][4k-2j]/9 

ee J[a{A) ■ cr{B)\ , (6) 

where 

J = (4k - 2j)/9 . (7) 

One sees that the effective exchange interaction be- 
tween two nearest neighboring trimers is antifer- 
romagnetic if 2k — j > and is ferromagnetic if 
2k — j < 0?i£ Thus the trimer-trimer interaction 



can be ferromagnetic even if all the spin-spin in- 
teractions are positive (antiferromagnetic) provid- 
ing j > 2k. These configurations are shown in Fig. 
[3J The elementary excitations within the ground 
manifold are those of a rectangular centered lattice. 
Then, if the trimers are antiferromagnetically cou- 
pled, standard spin-wave theory^ gives the doubly 
degenerate spin- wave energy w±(q) as a function 
of wave vector q, for — 7r/(2s) < q y < 7f/(2s) and 
-tt/(2\/3s) < q z < tt/(2\/3s), as 

w(q) = zJS^/T^W , (8) 

where z = 4 is the number of nearest neighbors, 
S =1/2, and 

7(q) = (lA)^cxp(iq- S) 
S 

= COs(sq y ) COs(v / 3sq 2 ) . (9) 

Here 8 is summed over nearest neighbor vectors be- 
tween trimers and s is the nearest neighbor sepa- 
ration in the Kagome lattice, as in Fig. 1. If the 
trimers are ferromagnetically coupled, then one has 
two nondegenerate modes whose energy is given by 

w±(q) = z|J|S[l± 7 (q)] ■ (10) 

Here (in Fig. [3J and below we give results for J = 1 
for the antiferromagnetic configuration of trimers 
with j = 0.15 and k = 0.2 and for the ferromag- 
netic configuration with j = 0.2 and k = 0.05. Note 
that transverse (H — ) modes of the antiferromagnetic 
configuration of trimers are doubly degenerate for 
all wave vectors. Also here and below note that the 
spectrum is always two fold degenerate for wave vec- 
tors on the face of the Brillouin zone [k y = w/(2s)] 
due to the Kramers-like degeneracy from the two- 
fold screw axis.[20j 

IV. EXCITON SPECTRUM 

Now we turn to the excitations out of the ground 
state manifold. 



A. Manifold at Energy J for the Antiferro 
Configuration 

Here we treat the case of antiferromagnetic cou- 
pling ( J > 0). The situation for this manifold 
is more complicated than that for the ground man- 
ifold. For the ground manifold we could develop 
degenerate perturbation theory for the manifold of 
2 N states of the system of N trimers in which each 



5 




FIG. 3: (Color online) As Fig. 1, the ferromagnetic 
(F) and antiferromagnetic (AF) arrangement of trimers, 
with spin orientation indicated by + or — . The inset 
graph shows the phase diagram of the trimer system in 
the j-k plane. The F-AF phase boundary based on the 
Neel state is at j = k and according to the trimer calcu- 
lation is at j = 2k. The latter calculation is more nearly 
correct when J is large compared to j or k, whereas the 
former is more accurate when J is not large compared 
to j and k. 



trimer independently occupies one of its two degen- 
erate ground states. The result was embodied in an 
effective Hamiltonian in which the interactions be- 
tween nearest neighboring trimers was given by Eq. 
(|7|) . For excitations near energy J we might consider 
the manifold of states in which one trimer occupies 
one of the excited states of Table I and all the other 
N — 1 trimers are distributed over the two degen- 
erate ground states. Strictly speaking this involves 
the solution to a many-body problem for the states 
of a spin excitation within the ground manifold and 
an exciton at excitation energy J or 3J/2. We will 
not treat this system with this degree of sophistica- 
tion. Instead, we will treat the manifold of excited 
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FIG. 4: (Color online) Spectrum of excitation energy, 
w(q) within the ground manifold for wave vectors in spe- 
cial directions. Here and below we plot the spectra for 
J = 1 for wavevectors on the lines joining F and X, X 
and S, and S and T, where F = (0,0), X = [tt/(2s),0), 
and S = [7r/(2s),7r/(2- v / 3s). For (j = 0.15, k = 0.2) 
one has a ferro configuration of trimers and for (j = 0.2, 
k — 0.05) one has an antiferro configuration of trimers. 
The modes shown here appear only in the transverse 
(H — ) response function. Modes are nondegenerate un- 
less labeled "(2)" to indicate a two-fold degeneracy. 



states at relative energy J or 3 J/2 when all the back- 
ground trimers are confined to their broken symme- 
try ground state. Thus our treatment is limited to 
the range of temperature T for which kT <C 3 ■ We 
therefore introduce operators a^(R) which take the 
trimer at R from its ground state to its nth excited 
state, where the labeling of sites is given in the first 
column of Table |U The Hamiltonian which describes 
the manifold at energy J is 



(ii) 



R 71=1 



where R is summed over trimer sites and n n = a^a n . 
Within the manifold near energy J the term in Eq. 

proportional to J is a constant and the nature 
of the band states is determined solely by the per- 
turbation V(J), which contains only terms propor- 
tional to j or k. To obtain results to leading order in 
the expansion parameter x, the perturbation V(J) 
is thus restricted to terms which conserve the un- 
perturbed energy J. Accordingly, the most general 
such form of V(J) is 



v(J) = £ 



2 

E - 

R, : R/ n.m— 1 



l (R,R')4(R)a m (R) 4(12) 



where the dots denote terms containing p creation 
operators (all with indices in the range 1,2) and p 
analogous destruction operators. Since we only con- 
sider nearest-neighbor interactions, we set 



R' — R n — R + d n 



(13) 



where 5 n are the nearest neighbor intertrimer dis- 
placements shown in Fig. [TJ 

<$i = — sj — V3sk , <5 2 = sj — V3sk 

S 3 = sj + V3sk, d 4 = -sj + \Z~3sk . (14) 

The effect of these 2pth order terms in Eq. (TT21 
on the mode energies is proportional to the (p — l)th 
power of the density of excitations. Since we as- 
sume that kT <C J, this density is small and we 
keep only the terms with p = 1. In addition we ig- 
nore the kinematic constraint which allows one to 
map the finite number of trimer states onto the in- 
finite number of bosonic states. [2l[ The discussion 
for the band at energy 3 J/ 2 is completely analogous 
to that for energy J and the analogous result holds 
for that case. So the band states are completely de- 
termined by the matrix c n!m (R, R'), or, as will shall 
see, by its Fourier transform which is a 4x4 matrix 
for the band at energy J and an 8 x 8 matrix for the 
band at energy 3J/2. To explicitly determine V{J) 
we must express the spin Hamiltonian in terms of 
the creation and annihilation operators of Eq. (|T2|) . 
The spin interaction between the fcth spin of an up 
trimer at R and the fc'th spin of a down trimer at 
R is 

S(fc; R) • S(fc'; R) = S z (k; R)S z (k'; R) 
+ [5+(fc; R)5_(fc'; R) + S_(fc; R)S + (k'; R)P25) 

Since 5+ and S- each involve at least one creation or 
annihilation operator, to construct the boson Hamil- 
tonian, we need only keep terms in these operators 
which are linear in the creation or destruction oper- 
ators. In contrast, since S z has a nonzero value in 
the ground state, we also need to keep terms in S z 
which involve one creation operator and one destruc- 
tion operator within the band. These considerations 
will be used implicitly below to limit the complexity 
of the mapping from spins to bosons. 

For the case of an "up" trimer at R (one whose 
ground state has S z = 1/2 and which we refer to as 
an "A" trimer) we find (keeping only terms linear in 
the boson operators) that 

S_(a;R) = at(R)/V3 - a|(R)/\/l8 

+a 6 (R)/V6 
S-(b;TL) = 24(R)/\/l8-2a 6 (R)/V6 
S_(c;R) = -a t 1 (R)/V3-4(R)/v / 18 

+a 6 (R)/\/6 . (16) 

The expression for S+(fe, R) are obtained by Hermi- 
tian conjugation. To determine the bosonic equiva- 
lent of S z we write 

S z = a + J2nm a nm4 a "» . (17) 



To determine the coefficients we require that the two 
representations lead to the same matrix elements. 
Thus if labels the ground state (i. e. whichever 
of the — J states of the trimer is the ground state), 
then, by taking matrix elements of both sides of Eq. 
(fT7|) we get 

a<3 = (0\S Z \0) , a n>m = (n\S z \m) , n^m 
a + a„. n = (n\S z \n) , n^O . (18) 

So for diagonal elements we must remember to sub- 
tract off the ground state value when identifying the 
bosonic matrix elements a nn . Thus 

S z (a;R) = l/3 + 4(R)/v / 12-4(R)/\/T8 
+a 2 (R)/\/l2- a 5 (R)/v / 18 - [2ni(R) + 2n 2 (R) 
+5n 3 (R) + 3n 4 (R) + n B (R) - ne(R)]/6 , 
S z (b; R) = -l/6 + 24(R)/Vl8 

+2a 5 (R)/v / l8- ni(R)/3 + 2n 2 (R) - n 3 (R)/3 
+n 5 (R)/3 + 2n 6 (R)/3 , 
S»(qR) = l/3-4(R)/VT2-4(R)/Vl8 
-a 2 (R)/VT2 - a 5 (R)/ v / 18 - [2m(R) + 2n 2 (R) 
+5n 3 (R) + 3n 4 (R) + n 5 (R) - n 6 (R) . (19) 

Here we needed to keep a^a p = n p terms in view of 
Eqs. (US) and (p~8|) . 

For the case of a "down" trimer (one whose ground 
state has S z — —1/2 at R and which we refer to as 
a "B" trimer) we similarly find that 

S+(a;R)=a t 2 (R)/V3 + al(R)/VT8 
-a 3 (R)/V6 

S+(b: R) = -24(R)/Vl8 + 2a 3 (R)/V6 
5+(c;R) = -al{R)/V3 + 4(R)/Vl8 

-o 3 (R)/Ve , 

S z (a]R) = -1/3- a\(R)/Vl2 
-a\{R)/Vl8-ai(R)/VT2-a i (R)/Vl8 
+ni(R)/3 + n 2 (R)/3 - n 3 (R)/6 
+n 4 (R)/6 + n 5 (R)/2 + 5n 6 (R)/6 , 
S z (b;R) = 1/6 + 24(R)/n/18 
+2a 4 (R)/Vl8- 2ni(R)/3 + n 2 (R)/3 
-2n 3 (R)/3 - n 4 (R)/3 + n 6 (R)/3 , 
S z (c;R) = -l/3 + 4(R)/Vl2 
-a\(R)/Vm + ai (R)/VT2 - a 4 (R)/\/l8 
+ni(R)/3 + n 2 (R)/3 - n 3 (R)/6 
+n 4 (R)/6 + » 5 (R)/2 + 5n 6 (R)/6 . (20) 

The next step is to write the interaction between 
trimers in terms of boson operators. Since we treat 
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here the case when the trimers are antiferromag- 
netically coupled, all interactions couple an up (A) 
trimer to a down (B) trimer. Since we treat only 
nearest neighbor interactions, we need consider only 
interactions between an up trimer at R and one of 
its four down neighbors at R± 6% and R ± 82- For 
the excitations band near energy J the boson Hamil- 
tonian is obtained in Appendix A. We define the 
Fourier transformed variables as 



4 >A (k) = iV-V3£ e* R 4(R) 
ReA 

ai !B (k)=iV- 1 /2^ e *R a t (R); (21) 
ReB 

where N is the total number of unit cells in the sys- 
tem. The quadratic Hamiltonian is of the canonical 
form: % — ^ q % q , where q is the wave vector and 
because we need consider only terms which conserve 
the unperturbed energy J, 



(22) 



where s = (n, r) and t = (n', r'). 

According to Table U excitations near energy J 
involve states 1 and 2 of the two spins in the unit 
cell, whereas excitations near energy 3J/2 involve 
states 3, 4, 5, snf 6 of the two spins in the unit cell. 
For excitations near energy J we write 



A = Jl + fcA fe + jAj , 

where X is the 4x4 unit matrix and Eq. 
Appendix A implies that 



A fc 
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4 3 7 (q) 

3 7 (q) 4 
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(23) 
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(25) 



The rows and columns of the matrices A are labeled 
in the order (1,A), (2, A), (1,5), (2,B). 

Thus the creation operators for the normal modes 
are ai i A(q) t , a2,s(q) t , and 

pi = [ai, B (q)t± 02 , A (q)t]/ > /2 , (26) 

with associated eigenenergies 

wiA(q) = ^2s(q) = J + (4fc + j)/9 (27) 
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FIG. 5: (Color online) As Fig. [4] but for excitations in 
the manifold near energy J for the antiferro (left) and 
ferro (right) configurations. The curve labeled "H — " in- 
dicates the energy in the transverse (H — ) response func- 
tion and those labeled "zz" are the energies in the longi- 
tudinal (zz) response function. The numbers in paren- 
theses indicate the degeneracy of the mode. 



and 



w±(q) = J+(4fc-5i)/9±(fc/3) 7 (q) . (28) 



These results are shown in Fig. [S] For an A (up) 
trimer a\ A corresponds to S- and a) 2 g corresponds 
to S + for a B (down) trimer. So these operators 
create transverse excitations and similarly one sees 
that p± create a longitudinal excitation. It may be 
surprising that, unlike for a Neel antiferromagnet, 
the longitudinal excitations exhibit dispersion, but 
the transverse ones do not. However, note that for 
a Neel antiferromagnet the dispersion comes from 
terms which here are eliminated since they do 
not conserve the large unerperturbed energy. 



B. Manifold at Energy J for the Ferro 
Configuration 

The calculations for the ferro configuration (in 
which all trimers start in their 'up' ground state) 
are similar and are done in Appendix C. In terms of 
Fourier transformed variables Eq. (IC31I) implies, in 
the notation of Eq. (|23|) . that 
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where the rows and columns are labeled in the order 
(A, 1), (A, 2), (B, 1), (B,2). The eigenvalues give 
the mode energies: 



wi, 2 = ■/ — j/9 — 4fe/9 ± 2fc 7 (q)/3 
w 3i4 = J + 5j/9-4£;/9±/c 7 (q)/3 



(31) 



These mode energies are shown for high symmetry 
wave vectors in Fig. [5j Since a' Al or a Bi connects 
the up ground state to a state with S z — —1/2, these 
operators correspond to an Thus we identify 

wi,2 as energies of transverse excitations and W3.4 as 
energies of longitudinal excitations as indicated in 
Fig. El 



C. Manifold at Energy 3J/2 for the Antiferro 
Configuration 

Here we adopt the same simplified approximation 
in which trimers not in excited states remain in their 
Neel state. Then, to leading order in the inter-trimer 
interactions, we only keep terms which are quadratic 
in the variables 3, 4, 5, and 5 and which conserve the 
total number of excitations. Thus analagously to Eq. 
(fTTjl we write 

6 

H(3J/2) = (3J/2)^^n„(R) + V(3J/2)(32) 

R ra=3 

The evaluation of V(3J/2) for the antiferro configu- 
ration is given in Eq. (|B38|) of Appendix B. In the 
notation of Eq. (j2"3")l , where we label the rows and 
columns of the matrices in the order 3 A, 6B, 4 A, 3B, 
5A, 4B, 6A, 5B, that result implies that 



A k = - 
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FIG. 6: (Color online) As Fig. [5] but for modes near 
energy 3 J/2 for the antiferro configuration. The highest- 
energy mode is not accessible in linear (in S) response 
theory. 



where 



T = -4V3 7 (q) , U = -8 7 (q) 
X = V3 7 (q) , Y = 2 7 (q) . 



(35) 



Thus we have the mode energies, with their degen- 
eracies in parentheses: 

ui = 3J/2 + 10fe/9- j/18 (2) , 

uj 2:3 = 3J/2 + 2fc/9-5j7l8±(4j-2fc) 7 (q)/9 (1) 

uj 4}5 = 3J/2 + 2fc/9-5j/18 

± V(4£ + j) 2 +3(fc-2j) 2 7 (q) 2 /9 (2) . (36) 

We determine the polarization of the modes as fol- 
lows. The mode lj\ involves state 3A which has 
Sa, z = —3/2 or state 6B which has Sb, z — 3/2 and 
is therefore not accessible via a single spin operator 
from the Sa,z = 1/2, Sb,z = —1/2 ground state. 
The modes L02 and L03 arise from states 5 A and 4B. 
State 5A has Sa.z = 1/2, which is activated from 
the Sa.z — 1/2 ground state by an Sa, z operator 
and state 4B has Sb, z = —1/2 which is activated 
from the Sb.z — —1/2 ground state by an Sb.z op- 
erator. The modes W4 and W5 arise from states 4A, 
3B, 6A, or 5B. State 4A has S A , Z = -1/2, which 
is activated by an Sa, - operator and state 5B has 
Sb.z = 1/2 which is activated by an Sb, + operator. 
States 3B or 6A lead to similar results. These modes 
(with their polarizations) are shown in Fig. [6l 



D. Manifold at energy 3 J/2 for Ferro 
Configuration 

The result of the calculation of I^(3J/2) for the 
ferro configuration is given in Eq. (|D44[) of Ap- 
pendix D, which implies, in the notation of Eq. (|23[) . 
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(38) 



where the rows and column are labeled in the order 
3A, 4A, 5A, 6A, 3B, 4B, 5B, 6B. We thereby find 
the mode energies to be 



Wl,2 


= 1.5J + 


W3,4 


= 1.5J + 


^5,6 


= 1.5J + 


^7,8 


= 1.5 J + 



j - 20k 

18 ' 
3j - 12k k - 2j 

18 9 

5j -4k 2k- Aj 



± 

18 9 

7j + 4k k - 2j 

18 - ~ 



7(q) 
7(q) 



7(q) ■ (39) 



The determination of the polarization of the mode 
is done as we did for the modes of Eq. (1551) . The 
results are shown in Fig. [7J 



V. NEEL SPIN WAVES 

In this section we compare the results obtained 
above with those from ordinary spin-wave theory. 
In Fig. [3] we show the 6 branches of transverse exci- 
tations from the Neel ground state. 

Note that apart from the lowest manifold, the two 
approaches lead to quite different spectra. As we 
showed above, the lowest manifold of trimer excita- 
tions is obtained by an exact mapping onto a Neel 
spin spectrum. One sees that for the anti configu- 
ration the energy scale of the lowest branch of spin 
waves is significantly larger for the trimer approach 
than for the Neel approach. This is because the 
trimer approach takes better account of quantum 
zero point motion that does the Neel approach. It is 



FIG. 7: (Color online) As Fig. [5] but for excitations in 
the manifold near energy 1.5J for the ferro configuration 
from Eq. (|39|l . The lowest-energy mode is not accessible 
in linear (in S) response theory. 




FIG. 8: (Color online) As Fig. Q] Neel (transverse) 
spin-wave spectrum for the anti (left) and ferro (right) 
configurations. Note that the two-fold degeneracy of the 
antiferromagnetic spectrum is broken along the low sym- 
metry S-F line. 



well known that zero point fluctuations tend to in- 
crease the spin- wave energies. This is shown by ex- 
act calculations for one dimensional spin chains [22j 
and by perturbative calculations for three dimen- 
sional systems. [23[ In contrast, for the ferro configu- 
ration the opposite effect occurs because the energies 
are proportional to the spin magnitudes. 



VI. NEUTRON DIFFRACTION 

Some aspects of neutron diffraction have been dis- 
cussed by Furrer et al£ and by Qiu et alii Here we 
discuss briefly the difference between the diffraction 
spectrum of the trimer system and that of the as- 
sociated Neel state. The elastic magnetic scattering 
intensity is proportional to 



|Q-F(Q)| 2 )5(Q-G), (40) 
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where G is summed over all reciprocal lattice vectors 
and the magnetic vector structure factor F is 

F(Q) « £<S.)e^, (41) 

T 

where r are the copper spin-positions given in Ta- 
ble HI1 and (S r ) is the thermal average of the spin at 
site t. For the Neel model, we take the spin-values 
as 0.5 while for the trimer model it is 1/6 and 1/3 
as shown in Table |TTJ To simplify the presentation 
we do not discuss the atomic form factor and the 
Debye- Waller factor. The magnetic elastic diffrac- 
tion intensities (apart from the thermal and mag- 
netic form factors) are summarized in Figure [9] for 
different collinear spin configurations along the crys- 
tal axes for both the trimer and Neel models, includ- 
ing AF and Ferro spin configurations. As expected, 
there are significant differences between the anti- 
ferro and ferro ordered trimer configurations. Also, 
for a given spin-configuration, the trimer model is 
significantly different than the Neel model. Due to 
smaller spin values in the trimer phase, the intenties 
are much weaker. Hence, observation of the mag- 
netic Bragg peaks would be much more difficult in 
the trimer phase than for the Neel model. Other 
than this difference, there are other differences at 
various scattering angle and it may be possible to 
distinguish the Neel and trimer model experimen- 
tally. In Figure [9j we also show nuclear scattering, 
which has some overlap with the strongest magnetic 
peaks. The unique magnetic peaks are at low angle 
and there are only a few of them. 




Q(A"') Q(A"') 



FIG. 9: Elastic magnetic Bragg peaks for different spin 
alignments for trimer and Neel models for the AF (left) 
and Ferro (right) configurations, respectively. The gray 
lines in the background show the nuclear scattering. 



TABLE II: Structure parameters [2^] (very similar re- 
sults are given in Ref. l24h for the distorted Kagome 
system Cu^OD^C l 11 ' 12 . Here x, y, and z are the P2ic 
fractional coordinates with respect to axes a = 9.1056A, 
b = 6.8151A, and c = 11.829A, with /3 = 30.825°. We 
choose the P2ic setting because the Kagome plane is 
x ~ in the this setting. The last column, S, shows the 
non-zero spin component in the trimer phase and is taken 
as along the a—, b—, and c— axis respectively in Figur({9] 
For the Neel model, we set the spin magnitude to 0.5 in- 
stead of 1/3 and 1/6. The CU3 sites are in the triangular 
lattice planes which interleave the Kagome planes, but 
our calculations do not include their moments. 
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VII. INELASTIC SCATTERING 
SCATTERING CROSS SECTION 

In this section we evaluate the inelastic cross sec- 
tion for the anti configuration at zero temperature. 
To do this we will construct the appropriate response 
functions, namely 

{{A-B)) = ^2{0\A\ri){n\B\Q)S(E n - tiu) ,(42) 

n 

where |0) denotes the ground state and the sum is 
over all states |n) with energy E n relative to the 
ground state. Here the operators A and B are pro- 
portional to the Fourier transforms of the spin oper- 
ators. In particular we will need 

S a p(q,u,) = ((S Q (q);^(-q))) . (43) 

Thus 

S+_(q,u;) = ^riHMqJIO)! 2 ^-?^) , 

n 

S_+(q,u;) = ^|(n|S + (q)|0)| 2 <S(i?„-M , 

n 

S zz (q,io) = ^|H^(q)|0)| 2 <S(£: n -M .(44) 

n 

(We later set Ti = 1.) To analyze the single-magnon 
contributions to these quantities we need to relate 
the spin operators to the normal mode operators. 
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Note that when we evaluate Eq. (|4"4")l at zero tem- 
perature, only contributions to the operator ^(q) 
proportional to creation operators are nonzero. We 
will quote results for the transverse and longitudinal 
cross sections, given respectively by 

/trans (q, w) = <S+- (q, ui) + <S_+ (q, w) , 
7ion g (q,w) = S zz (q,uj) . (45) 

In the calculations which follow we use the notation 
introduced in Sec. VI. 



1. GROUND-STATE EXCITATIONS 

We first consider inelastic scattering from pseudo- 
spin waves. Accordingly, we discuss spin- wave the- 
ory for this situation. We express the pseudo-spin 
operators in terms of boson creation operators, 
and cjg for the A (up) and B (down) trimers, respec- 
tively, as 



a z (A) = 1/2 ~c\c A 
<j z {B) = -1/2 + c A c A 

and (with a± — <r x ± icr y ) 



CB 



(46) 



(47) 



Then, following the standard spin-wave treatment 
for such a spin 1/2 system we write 



c A (q) = N- 1 ' 2 ]T e-^ R c A (R) 



(48) 



Re a 



and similarly for c B (q), where R is summed over 
all the N positions of A trimers. Then the boson 
Hamiltonian H = X) q ^q a ^ quadratic order is 

= 2 l 7{c+ l (q)c^(q) + 4(-q)c B (-q) 
+7(q)[4(q)4(-q) + CA(q)c B (-q)]} .(49) 

Then the operators which create normal modes are 
pt(q) and rf(q), which are determined by 



cUq) = '(q)p f (q) -fra(qM-q) 



(50) 



and 



where 

^(q) 2 
Z(q)m(q) 



cs(-q) = -™(q)p f (q)+^(q)?7(-q) , (51) 

l-e(q) 



1 + e(q) , ^2 

m(q) 

7(q) 

2e(q) ' 



2e(q) 



(52) 



with e(q) = [1 — 7(q) 2 ] 1 ^ 2 - Apart from the constant 
zero point energy one has 

U =^a;(q)[ /9 t( q ) p (q) + r? t (q)7? ( q)] , (53) 
q 

where Eq. © is w(q) = 2.7e(q). 

Using Eq. (jj) we note that the Fourier transform 
of the spin operators is 

ReA 

+N- 1 ' 2 a tt (R)e- 4q - R r s (q) .(54) 
Res 



Here we have introduced the trimer form factors 

(55) 



rx(q) = -cos(q-nx) - g , 



where "fix incorporates the locations of the sites of 
trimer X relative to its center of gravity: 

n A = 5(0,1/2,73/2) 

h B = s(0,-l/2,V3/2) . (56) 

Within the approximation of a Neel state 

rx(q) = 2cos(q • h x ) - 1 • (57) 

When B in Eq. (|42)) is proportional to 5_(q) wc 
have (at zero temperature) 

B = r A (q)c^(q) +r B (q)c B (-q) 

-> [/(q)r A (q) - T B (q)m(q)]pt( q ) + . . . (58) 

where the dots indiocate terms involving r/(q) which 
do not contribute at zero temperature. In /trans we 
also have the contribution when B in Eq. (|42|) is 
proportional to S+(q), in which case 

B = [Kq)r B (q)-r A (q)m(q)]r ? t(- q ) + ...(59) 

Thus the contribution to the inelastic transverse 
cross section is given by 

Itrans(qw) = { [*(q) 2 + m(q) 2 } [r A (q) 2 + t b (q) 2 ] 

-4?(q)m(q)T j4 (q)T B (q)}5[w - w(q)] 

= ^y{Mq) 2 +T B (q) 2 ] 

-2 7 (q)r A (q)r B (q)}5[a; - w(q)](60) 

In the case of a standard two-sublattice antiferro- 
magnet, one has the same result but with 
T A(q) = r s(q) = 1. In that case inelastic scattering 
cross section for spin waves alternates in intensity as 
one goes from one Brillouin zone to the next due to 
the alternating sign of 7(q). Here the result is more 
complicated because of the form factor of the unit 
cell, reflected by the factor rx(q). 
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2. EXCITONS NEAR ENERGY J 

To get the response near energy J for the anti- 
ferro case we need to construct the nonzero matrix 
elements required to evaluate Eq. (|l4|l . To obtain 
the cross section near energy J we only consider con- 
tributions which involve a\ (R) or a\ (R) • From Eq. 
(|16p and following we see that the only nonzero con- 
tributions of this type are, 

S- (a, R) = 4(R)/V3 = -SL(c,R) 
S z (a,K) = 4(R)/VU= -S z (c,TL) , (61) 
where R is an A, or up, trimer and 

5+(a,R) = 4( R )/^3 = -5+(c,R) 
5 2 (a,R) = -4(R)/\/l2 = -S z (c,R) , (62) 

when R is a B, or down, trimer. These results lead 
to 

5_(q) - aUq)[e™- e - lq - n -1/\/3 

= (2i/N/3)a(q)4 A (q) , (63) 
where £x(q) = sin(q • n^). Similarly 
5+(q) - (2»/V3]e B (q)4 fl (q) , 

s,(q) = ( i /V3)[a(q)4A(q)-^(q)4s(q)I64) 

Then, using Eq. (|44l) . we have 

/trans = (4& (q) 2 /3) ( (a 1A ; 4a) ) 

+ (46 3 (q) 2 /3)((a 2B ;4 i3 )) , (65) 

where Eq. (|2~7| gives 

((aiA;a| A )) = ((a2s;a^ B )) = <5[cj - cjiA(q)](66) 

Also Eq. (gll) gives 

4s (q) = [pUq) + pL(q) /v 7 ^ 

a + 2A (q) = ^(q)-pl(q) /V5 , (67) 

so that 

5*(q) = (*/V6)^eA(q)-efl(q)]pUq) 

-[6t(q)+6*(q)]/?l(q)]V (68) 

Then we obtain 

«5,(q);S.(-q)» 

= (l/2)[^(q)+eB(q)] 2 ((p + (q);pUq)» 

+ (l/2)[^(q)-^(q)] 2 ((p_(q);p t _(q))) ,(69) 

where Eq. (|2"5)l gives that 

((P ± (q);p f ±(q)» = *[w-w±(q)] . (70) 



3. EXCITONS NEAR ENERGY 3 J/2 

Here we keep only contributions involving creation 
operators <r n , with n > 2. In this case 

5_(a,R) = -4a/V18 = 5- (c,R) 
5_(6,R) = 24 A / v / l8 
5+(a,R) - 4 A /V6 = 5+(c,R) 
5+(6,R) = -24 A /V6 

5 z (a,R) = -«4 A /Vl8 = £L(c,R) 

5,(6, R) = 24^/\/l8 

(71) 

when R is an A, up, site. Also 

5+(a,R) = a\ B /Vw = 5_(c,R) 
5+(6,R) = -2a\ B /V\% 
5_(a,R) = -4 B /V6 = 5_(c,R) 
5_(6,R) = 2al B /V6 
S z (a,R) = -4 B /v / 18 = 5_(c,R) 
S z (b,R) = 2a\ B /VT8 

(72) 

when R is an B, down, site. Thus 

5_(q) = -M A (q)4 A (q)/Vl8 
-/i B (q)4 s (q)/v / 6 

5+(q) = MA(q)4yi(q)/^ / 6 

-/i B (q)4 B (q)/Vi8 

5,(q) = -[MA(q)4A(q) + ^(q)4 B (q)]/v^3) 

where 

Mx(q) = 2cos(q- n x ) - 1 . (74) 

The intensities can be obtained by inverting the 
transformation which diagonalizes V(3J/2) whose 
eigenvalues are given in Eq. (j36]). Since the al- 
gebraic expression for the mode intensities are too 
complicated to be enlightening, we confine ourselves 
to some general remarks. We verify that 5_(q) in- 
volves the third and fourth rows and columns of the 
dynamical matrices of Eqs. (|33p and (pHf . Like- 
wise 5+ (q) involves the seventh and eighth rows and 
columns of the dynamical matrices of Eqs. (|33[) and 
Thus the transverse response is associated with 
modes u>4, and ws, in agreement with our previous 
identification. Similarly we confirm the identifica- 
tion of u>2 and uj^ as belonging to the longitudinal 
response. 
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VIII. CONCLUSIONS 

We may summarize our results as follows. 

A) The lowest energy modes of the trimer system 
shown in Fig. |4]arc only slightly different from what 
one gets (see Fig. [5J using the Neel approximation 
for the ground state. There is a slight difference in 
symmetry in that the breaking of degeneracy of Neel 
spin wave in nonspecial directions does not occur 
in leading order of perturbation theory within the 
trimer approximation. 

B) The elastic diffraction pattern shows differences 
(see Fig. ^ which, in principle, allow one to distin- 
guish between a trimer system and one that is closer 
to the Neel limit. 

C) The excitation spectra at high energy we have ob- 
tained show dramatic differences between the trimer 
and Neel limits. In the former case, well defined 
modes appear in the longitudinal response functions. 
In general the trimer limit gives rise to many more el- 
ementary excitations and thereby provides a conclu- 
sive way to identify a system as being in the trimer 
limit. 

D) A possible future project would be to develop 
an interpolation scheme to pass between the quali- 
tatively different Neel and trimer limits. 

ACKNOWLEDGMENTS. ABH was supported in 
part by a grant from the department of commerce. 



Thus this interaction leads to the Hamiltonian 

U = j[-2m(S 3 )+n 2 (8 a )}/9 

-j[ni(0)+na(0)]/18. (A5) 

2. a at interacts with c at 83 and 84 

Here 

S-(a) = a\(0)/V3, 5+(a) = ai(0)/V3, (A6) 
a / \ 1 4(0) a 2 (0) m(0) ng(Ql 

and, where 5 assumes the values #3 and S4, 

S-(c,5) = -a 2 {d)/V3, (A8) 
S+(c,6) = -<4(<S)/V3, (A9) 
S z (c,6) = -l/3 + al(S)/Vu + ai(S)/Vl2 

+n 1 (S)/3 + n 2 (6)/3 . (A10) 

Thus this interaction leads to the Hamiltonian 

n = §5^M*) + na(*) + ni(0) + na(0)] 
S 

+^^2[4(0)ax(8) + al(8)a 2 (0)} .(All) 



Appendix A: Antiferro Excitations at Energy J 



3. b at interacts with a at £2 



Here position coordinates are given relative to R 
a lattice site occupied by an 'up' trimer. Thus a 2 (0) 
denotes a2(R), ai(^i) denotes ai(R + 5i) and so 
forth. We treat the interaction of one of the spins 
(a, b, or c) of the trimer at R with one of the spins 
(a, b, or c) of a neighboring trimer at R+<5„, for n 
1,2,3,4. In this section we drop all terms referring 
to states n > 2 since such states occur at energy 
3J/2. Also, as mentioned, we drop all terms which 
are off diagonal in J. 



1. a at interacts with b at #3 

Within the band at energy J we may write 

S-(a) = oJ(0)/V3, 5 + (o) = oi(0)/v'3,(Al) 

S z (a) = I + 4W + ^_M0)_^ 2) 
3 y/12 VT2 3 3 V ' 



and 



S-(b,S 3 ) = 



S+(M 3 ) = (A3) 



S.(Ms) = i-^*2) + ^. (A4) 

6 3 3 



Here 



S- (b) =0 S+{b)=0 (A12) 
S z (b) = -l/6 + 2n 2 (0)/3-m(0)/3 (A13) 



— and 



S-(a, 8 2 



V3 ' V3 
S z (a,S 2 ) = -1/3 - a\(S 2 )/Vl2 - ai(S 2 )/Vu 
+n 1 (S 2 )/3 + n 2 (S 2 )/3 . (A15) 

These interactions lead to the Hamiltonian 

H = -j[n 2 (6 2 ) + ni (6 2 )}/18 

+j[ni(0)-2n 2 (0)]/9. (A16) 

4. b at interacts with c at #4 



Here 



S-(b) = , S+(b) = , (A17) 

S z (b) = -l/6 + 2n a (0)/3-n 1 (0)/3 (A18) 
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and 
S-(c,S 4 ) 



a2{Si) S + (c,S 4 ) = -^Am 



V3 ' ~^~""» " V3 V 
S z (c,5 4 ) = -1/3 + a\(S 4 )/VT2 + ai (S 4 )/VT2 
+n 1 (d 4 )/3 + n 2 (d 4 )/3 . (A20) 

These interactions lead to the Hamiltonian 

H = -j[ni(5 4 )+n 2 (« 4 )]/18 

+j[m(0)-2n 2 (0)]/9 . (A21) 

5. c at interacts with a at <5i and 82 

Here 

S_(c) = -a\{0)/\/3, S+(c) = -oi(0)/>/3(A22) 

s, ( c) = I-^-^-^M-^) 

U 3 712 712 3 3 V ; 



and, where <5 assumes the values <5i and <5 2 , 

S-(a,S) =0, S+(6,5) = 0, (A24) 
S z (a,S) = -l/3-4(5)/\/l2-ai((5)/\/l2 

+m(5)/3 + n 2 ((5)/3 . (A25) 

These interactions lead to the Hamiltonian 

U = \ EM'S) + n 2 (d) + m(0) + n 2 (0)] 
9 5 

+ ^£[4(0KW+4(<5)a 2 (0)]. (A26) 



7. Summary 

Summing all the above contributions we get the 
Hamiltonian for the band at energy J for the anti- 
fcrro configuration 

H(J) = ^(j[n 1 (R)-5n 2 (R)+n 2 (R 1 ) 

R \ 

-5m(Ri)]/9 + 4fc[m(R) + n 2 (R) 
+m(Ri) + n 2 (Ri)]/9 
+ ^A;[4(R)ai(R + <5) 



+4(R + <5)a 2 (R)]/12 



(A31) 



where S is summed over the four values shown in 
Fig. 1. 



Appendix B: Antiferro Excitations at Energy 

3J/2 

1. a at interacts with b at 83 



Here 



S-(a) = -4(0)/\/l8 + a 6 (0)/V6 , 
S+(a) = -a 4 (0)/\/T8 + 4(0)/v / 6, (Bl) 
1 4(0) 05(0) , ne(0) 



S z (a) 



3 ^18 Vl8 6 
n 5 (0) n 4 (0) 5n 3 (0) 



• (B2) 



6. c at interacts with b at <5i 



and 



Here 



S-(c) = -4(0)/V3, S+(c) = -ai(0)/v^A27) 
^ (c) = 1 4(0) o 2 (0) , ni(0) , n 2 (0 A ) n 



and 



3 JU Vl2 3 



SI (Mi) - 0, S+(Mi) = 0, 



+ + ^|A28) 



1 2m(gi) n 2 (Ji) 
S*(Mi) = + — - — . (A29) 

These lead to the Hamiltonian 

H = -j[m(0) + n 2 (0)]/18 

+j[-2m(5i)+n 2 (*i)]/9 . (A30) 



S-(b,S 3 ) = -2a 5 (S 3 )/Vl8 + 2al(S 3 )/V6,(B3) 
S+(b,S 3 ) = -2al(S 3 )/Vl8 + 2a 3 (S 3 )/V6,(B4) 

s z ( b ,s 3 ) - 1 ■ 2fl ^) ■ 2a4(<53) 



6 vT8 V18 
Mh) _ n 4 (S 3 ) _ 2n 3 (8 3 ) 
3 3 3 k ' 

These interactions give rise to the Hamiltonian 

H = V3j[-a\(0)a 3 (S 3 ) - 4(8 3 )a 6 (0) 
-4(5 3 )a 4 (0)-4(0)a 5 (5 3 )]/18 
+j[4n 6 (S 3 ) - 4n 4 (<5 3 ) - 8n 3 (<5 3 ) + n 6 (0) 
-n 5 (0) - 3n 4 (0) - 5n 3 (0)]/36 
+j[-4(0)a 4 ((5 3 ) - a 5 (0)4((5 3 )]/9 . (B6) 
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a at interacts with c at #3 and £4 



Thus the Hamiltonian for this interaction is 



Here 



S-(a) 
S+(a) 



-4(0)/VT8 + a 6 (0)/\/6 , 
-0,4 (0)/\/l8 + 4(0)/\/6 , 
4(0) a B (0) ne(0) 



^ VT8 718 ' 6 
n 5 (0) n 4 (0) 5n 3 (0) 
6 2 6 

and, with J = £3 or d = S 4 , we have 

S-(c,6) = a 5 (5)/VT8-4(5)/V6 , 
5+(c,<5) = 4W/^18- a 3(«5)/V6 , 



(B7) 
(B8) 



(B9) 

(BIO) 
(Bll) 



V ; 3 VI 



04(5) 5ng(5) 



| n 5 (g) n 4 (5) 



18 6 
n 3 (8) 



2 6 6 

These interactions lead to the Hamiltonian 



(B12) 



5 



^(a\(0)a 3 (5)+al(5)a 6 (0) 



+4(0)a 5 ((5) + 4(<*)a 4 (0) 



18 



5n 6 (<5) + 3n 5 (<5) + n 4 ((5) - n 3 (S) 



-n 6 (0) + n 5 (0) + 3n 4 (0) + 5n 3 (0) 



(B13) 



Here 



3. b at interacts with a at £2 

S-(b) = 24(0)/\/l8-2a 6 (0)/V6 , (B14) 
S+(b) = 2a 4 (0)/Vl8- 24(0)/\/6 , (B15) 

6 VT8 Vl8 



2n 6 (0) n 5 (0) n 3 (0) 
+ — 5 1 5 3~ l Bib ) 



and, 



S-(a,5 2 ) = a 5 (5 2 )/Vm-al(5 2 )/V6, 
S+(a,S 2 ) = al(S 2 )/VT8-a 3 (d 2 )/V6, (B17) 
1 4(^2) 04(^2) . 5n 6 (<5 2 ) 



5,(o,5 2 ) = -- ' 1 



+ 



18 V18 6 
n 5 (S 2 ) n 4 (S 2 ) n 3 (8 2 ) 



.(B18) 



ft = \/3j[-4(0)a 3 (<5 2 ) - 4(5 2 )a 6 (0) 
-4((5 2 )a 4 (0) - 4(0)a 5 (5 2 )]/18 
+j[-5n 6 (5 2 ) - 3n 5 ((5 2 ) - n 4 (5 2 ) + "3(^2) 
-8n 6 (0)-4n 5 (0)+4n 3 (0)]/36 
+j[-oj(0)a 4 (5 2 ) - 4(^)a 5 (0)]/9 . (B19) 



4. b at interacts with c at <?4 



Here; 



S_(6) = 24(0)/v/l8-2o 6 (0)/\/6 , (B20) 
S+(b) = 2a 4 (0)/\/l8- 24(0)/\/6 , (B21) 
1 24(0) 2 Q5 (0) 2n 6 (0) 

" "6 + T7iT + T7iT + ^- 



+ 



n 5 (0) n 3 (0) 



(B22) 



and 



S_(c,5 4 ) = a 5 (5 4 )/\/l8-4(<5 4 )/\/6, 
S+(c,6 4 ) = al{5 4 )/Vl8-a 3 (6 4 )/V6, 



S z (c,Si) 



1 _ 4(^4) _ ^4(^4) 5n 6 (d 4 ) 
3 ~ ^18 y/18 6 
+ n^4) + rc 4 (<S 4 ) _ n 3 (S 4 ) 



2 6 6 

These results lead to the Hamiltonian 

U = V3j[-a\(0)a 3 (S 4 ) - 4(<* 4 )a 6 (0) 
-4(5 4 )o4(0) - 4(0)a 5 (5 4 )]/18 
+j[-5n 6 (5 4 ) - 3n 5 (<5 4 ) - n 4 (5 4 ) + n 3 (<5 4 ) 
-8n 6 (0)-4n 5 (0)+4n 3 (0)]/36 
+j[-4(0)a 4 ((5 4 ) - 4(^)a 5 (0)]/9 . (B24) 

5. c at interacts with a at Si and 82 

Here 

S-(c) = -4(0)/Vl8 + a 6 (0)/v / 6 , (B25) 
S+(c) = -a 4 (0)/\/l8 + 4(0)/\/6 (B26) 
S Z ( C ) = 1 4(0} gg(0) , "e(0) 



3 Vl8 Vl8 6 
n 5 (0) n 4 (0) 5n 3 (0) 
6 2 6 

and, where <5 assumes the values Si and <5 2 , 
5_(a,5) = a 5 (d)/VT&- al(5)/V~6, 



(B27) 
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S+(a, S) = a 5 (S)/Vl8 — as(S)/ VQ, (B28) antiferro configuration 



5n 6 (d) n 5 (S) n 4 (S) _ n 3 (S) 
6 2 6 6 ' 

Thus the Hamiltonian from this interaction is 



(B29) 



s 



^(a\(0)a 3 (5)+ae(0)ai(8) 



+04(0)4(5) + aJ(0)a B («) 
fc 



18 



5n 6 (<5) + 3n 5 (5) + n 4 ((5) - n 3 (<5) 



-2n 6 (0) + 2n 5 (0) + 6n 4 (0) + 10n 3 (0) 



+4(0)a 4 ((5) + a 5 (0)4(5) 



(B30) 



Here 



6. c at interacts with b at <5i 

S*_(c) = -4(0)/\/l8 + a 6 (0)/V6 (B31) 
S+(c) = -a 4 (0)/\/i8 + 4(0)/\/6 (B32) 
a ( \ 1 4(0) o 6 (0) n 6 (0) 

^z(c) = 7=- + 



3 718 6 
n 5 (0) n 4 (0) 5n 3 (0) 



(B33) 



and 



S-(b,8!) = -2o 5 (5i)/n/18 + 24(«5i)/v^B34) 
S+(Mi) = -2aJ(«i)/v^8 + 2a 3 (5i)/v^B35) 
<? rfc -S ^ 1 4. 2fl j(^) , g04(gl) 



i »e(gi) fu(gi) 2n 3 (tfi) 

+ ~ 3 S - ' (B36) 

Thus the Hamiltonian from this interaction is 

ft = V3j[-a\{0)a 3 (S) - 4(S)a e (0) 
-4((5)a 4 (0)-4(0)a 5 (5)]/18 
+j[4n 6 (5) - 4n 4 (<5) - 8n 3 (<5) + n 6 (0) 
-n 5 (0) - 3n 4 (0) - 5n 3 (0)]/36 
+j[-4(0)a 4 (5) - 4(5)a 5 (0)]/9 . (B37) 

7. Summary 

Summing all the above contributions we get the 
Hamiltonian for the band at energy 3 J/2 for the 



R 



[-n 3 (R) - 3n 4 (R) - 5n 5 (R) - 7n 6 (R) 



-n 6 (Ri) - 3n 5 (Ri) - 5n 4 (R 4 ) - 7n 3 (Ri)]/18 
+ ( -[4( R )«4(R + 6) + 4(R + <5)a 5 (R)]/9 

s V 

-\73[4(R)as(R + 5) + 4(R + <5)a 4 (R) 



+4(R + S)a e (R) + 4(R)a 5 (R + S)]/I8j 



+*E 



R L 



2[-n 6 (R) + n B (R) + 3n 4 (R) + 5n 3 (R) 



+5n 6 (Rx) + 3n 5 (Ri) + n 4 (R 4 ) - n 3 (Ri)]/9 
+ E ( >/3[aJ(R)a 3 (R + 5) + J(R + <5)a 4 (R) 

+4(R + S)a e (R) + 4(R)a 5 (R + *)]/36 
+ [4(R)a 4 (R + S) + 4(R + «5)a 5 (R)]/18j438) 



Appendix C: Ferro Excitations at Energy J 
1. a at interacts with b at 83 

Here 

S-(a) = a\(0)/V3,S+(a) =ai(0)/V3, (CI) 
« , x 1 4(0) , 02(0) m(0) ngjO)^ 

Sz{a) 3 + 7lf + 7H~~ s* 3 2) 



and 



S.(b,S 3 ) = 0, S+(M 3 ) = 0, (C3) 
1 m(5 3 ) 2n 2 (<5 3 ) 



5,(6, 5 3 ) 



+ 



6 3 3 
These results lead to the Hamiltonian 



• (C4) 



ft = j[2n 2 (5 3 )-ni(« 3 )]/9 

+j[m(0)+n 2 (0)]/18 . (C5) 

2. a at interacts with c at 63 and <5 4 



Here 



S_(a) = 4(0)/V3, S+(o) = oi(0)/>/3, (C6) 
q (n , 1 4(0) , 02(0) m(0) n 2 (QL 
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and, where S assumes the values £3 and £4, 

(n r\ 4(g) o / ^ , P81 

6_(c,d) = — -j=-, b+{c,5) = — -^=- ,(C8) 

S z (c,«) = 1/3-4(5)/v^2-o 2 (5)/>/12 

-ni(5)/3-n 2 ((5)/3 . (C9) 

Thus, we have that 



k 



--(n 1 (S)+n 2 (S)+n 1 (0)+n 2 (0) 



+^{-4(S)a 2 (0)-4(S)a 2 (0) 



--( 4(0)01(5) + 4(5)ai(0) 



3. b at interacts with a at £2 



(CIO) 



Here 



S_(6) = , S+(b) = , 



(Cll) 



^) = -i + ^--M, (C12) 
J 3 



and 



6_(a,d 2 ) = -^= , fr + (a,d 2 ) = ^ (C13) 

S z (a,6 2 ) = l/3 + 4{S 2 )/VT2 + a 2 (S 2 )/VT2 

-ni((5 2 )/3-n 2 ((5 2 )/3 . (C14) 

Thus, we obtain the Hamiltonian 

^=jM<5 2 )+n 2 (5 2 )]/18 
+j[-m(0)+2n 2 (0)]/9 . (C15) 



Here 



4. b at interacts with c at 84 



S-(b) = , S+(b) = 0, (C16) 
+ (C17) 



and 



6-(c,d 4 ) = -^-,6+(c,d 4 ) = ^-(C18) 

S»(c,« 4 ) = l/3-a 2 (5 4 )/\/T2-a 2 (5 4 )/\/T2 

-m(5 4 )/3-n 2 (<5 4 )/3 . (C19) 

These interactions give rise to the Hamltonian 

H=j[n 1 (6 i )+n 2 (5 4 )}/18 
+j[-m(0) + 2n 2 (0)] 9 . (C20) 



5. c at interacts with a at Si and £2 

Here 



S-(c) = 



al(0) S + (c) = , (C21) 



V3 ' ^ ' V3 
S z (c) = l/3-4(0)/VT2-a 2 (0)/\/l2 

-m(0)/3-n 2 (0)/3 . (C22) 

and, where 5 assumes the values Si and S 2 , 

S-(a,S) = S +M = ^(C23) 

S z (a,<5) = l/3 + 4(5)/v / 12 + a 2 (5)/\/T2 

-m((5)/3-n 2 (5)/3 . (C24) 

These results lead to the Hamiltonian 



s 



^[n 1 (d)+n 2 (d)+n 1 (0)+n 2 (0) 



-^(4(0)a 2 (5) + 4Wa 2 (0) 
-^4(0)ai(5) + 4(5)ai(0)) - (C2r>) 



6. c at interacts with b at £1 



Here 

S-(c) - -^P, S + (c) = -^, (C26) 



V3 



V3 ' 



Q , x _ 1 4(0) 02(0) m(0) r^m 
Sz[c) 3"7Tf ~lm ~~ 3^ 7) 



and 

SI (Mi) = 0, S+(6,5i)=0 (C28) 

S,(Mi) = 1 + 2^*0 _ 21^1) (C29) 
3 3 

These terms give rise to the Hamiltonian 

H=j[m(0)+n 2 (0)]/18 
+j[2n 2 (5 1 )-n 1 (5 1 )]/9 . (C30) 



7. Summary 

Summing all the above contributions we get the 
Hamiltonian for the band at energy J for the ferro- 
conhguration as 



R 



j[5n 2 (R) - rci(R) - m(Ri) + 5n 2 (Ri) 



18 



-4fc[m(R) + n 2 (R) + m(Ri) + n 2 (Ri)]/9 
~J2 \ [2a t 1 (R)a 1 (R + S) + 4(R)a 2 (R + «J) 

12 s V 

+2a{(R + 5)ai(R) + 4( R + $)a 2 (R)^031) 



Appendix D: Ferro Excitations at Energy 3J/2 
1. a at interacts with b at S3 



Here 



S-(a) 
S+(a) 

S z (a) 



-4(0)/\/l8 + a 6 (0)/V6 , (Dl) 
-a 4 (0)/\/T8 + 4(0)/V6 , (D2) 
1 4(0) a 5 (0) ne(0) 



(D3) 



3 V18 V18 6 
n 5 (0) n 4 (0) 5n 3 (0) 



and 



S-(b,S 3 ) = 24(<5 3 )/\/l8-2a 6 ((53)/\/6 , (D4) 
S+(b,S 3 ) = 2a 4 ((5 3 )/\/T8-24((53)/\/6 , (D5) 



S,(M) = -l/6 + 2aJ(5 3 )/V^8 + 2a B (53)/>/l8 
+2n 6 (<5 3 )/3 + n 5 (<5 3 )/3 - n 3 (<5 3 )/3 .(D6) 

Thus these interactions lead to the Hamiltonian 

H = j[-al(0)a 4 (6 3 ) - 3oJ(5 3 )a«(0) 
-4(5 3 )a 4 (0) - 3<4(0)o6(5 3 )]/18 
+j[-4n 3 (5 3 ) + 4n 5 (<5 3 ) + 8n 6 (<5 3 ) 
-n 6 (0) + n 5 (0) + 3n 4 (0) + 5n 3 (0)]/36 
+i[-4(0)a 5 (5 3 ) - 4((5 3 )a 5 (0)]/9 . (D7) 

2. a at interacts with c at 63 and ($4 

Here 



S*_(«) = -4(0)/\/l8 + a 6 (0)/v / 6 , 
S+(a) = -a 4 (0)/\/l8 + 4(0)/\/6 , 
C , x 1 4(0) o 6 (0) n 6 (0) 

M } _ 3 17m Trf + 

n 6 (0) n 4 (0) 5n 3 (0) 



6 



6 



(D8) 
(D9) 



(D10) 



and, where 5 assumes the values 63 and 5 4 , 

S-(c,S) = -a\{6)/VT$ + a 6 (5)/V6 , (Dll) 



S+(c,5) = -a 4 (5)/\/T8 + 4(<5)/\/6 , (D12) 

Sz{c , 5) = 1 4g_^_^ D13) 
V ; 3 Vl8 >/l8 6 v ' 



m{8) n 6 {8) | n 6 (<5) 



Thus 



36 



oi(«)a 4 (0) + 3oj(5)a6(0) 



+4(0)a 4 (^)+34(0K(<5)) 



18 



-5n 3 (5) - 3n 4 (8) - n 5 (S) + n e (8) 



+n 6 (0) - n 5 (0) - 3n 4 (0) - 5n 3 (0) 



-^(4(0KW + 4(«53)a 5 (0)) . (D15) 



3. b at interacts with a at S 2 



S-(b) = 24(0)/VT8- 2a 6 (0)/V6 , (D16) 
S+(b) = 2a 4 (0)/Vl8- 2al(0)VE , (D17) 
S*(6) = -l/6 + 24(0)/\/l8 + 2a 5 (0)/\/T8 

+2n 6 (0)/3 + n 5 (0)/3 - n 3 (0)/3 (D18) 



and 



S_(M 2 ) = -4(5 2 )/v / T8 + a 6 (<5 2 )/\/6 , (D19) 
S+(a,<* 2 ) = -a 4 (d 2 )/VTs + al(S 2 )/V6 , (D20) 
1 _ 4(^2) _ 05(^2) _ 5n 3 ((5 2 ) 
3 v/T8 \/l8 6 

(D21) 



5 z (a,5 2 ) 



^4(^2) _ "-5(^2) ^5(^2) 
2 6 6' 

Thus, we obtain the Hamiltonian 

H = j[-at(0)a 4 (S 2 ) - 3at(S 2 )a 6 (0) 



-a\(d 2 )a 4 (0) - 3al(0)a 6 (8 2 ) 



/18 



+j[5n 3 (£ 2 ) + 3n 4 (<5 2 ) + n 5 (S 2 ) 
-n 6 (6 2 ) + 8n 6 (0) + 4n 5 (0) - 4n 3 (0)]/36 
+j[-4(0)a 5 (S 2 ) - 4(<5 2 )a 5 (0)]/9 . (D22) 

4. b at interacts with c at 64 



Here; 

S-(b) = 24(0)/\/l8- 2a 6 (0)/V6 



(D23) 
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S+(b) = 2a 4 (0)/\/l8- 24(0)/\/6 (D24) 
S z (b) = -l/6 + 24(0)/yi8 + 2a 5 (0)/\/l8 

+2n 6 (0)/3 + rc 5 (0)/3 - n 3 (0)/3 (D25) 



and 



S-(c,5 4 ) = -4(5 4 )/\/l8 + a 6 (5 4 )/\/6 , (D26) 
S+{c,S 4 ) = -a 4 (6 4 )/Vl8 + al(6 4 )/V6 , (D27) 



^ (C ' 54) = 3" V18 



1 4(5 4 ) a 5 (<5 4 ) 5n 3 (<5 4 ) 



18 6 

n 4 (5 4 ) n 5 (5 4 ) n 6 (<5 4 ) 



(D28) 



2 6 6 

Thus we obtain the Hamiltonian 

ft = j[-4(0)a 4 ((5 4 ) - 3aS(« 4 )o6(0) 
-a|(5 4 )a 4 (0) - 34(0)a 6 ((5 4 )]/18 
+j[5n 3 (<5 4 ) + 3n 4 (<5 4 ) + rc 5 ((5 4 ) - n 6 (<5 4 ) 
+8n 6 (0) + 4n 5 (0) - 4n 3 (0)]/36 
+j[-4(0)a 5 (5 4 ) - 4(5 4 )a 5 (0)]/9 . (D29) 

5. c at interacts with a at <5i and S2 



Here 

S-(c) 

S+(c) 

Sz(c) 



-4(0)/Vl8 + a 6 (0)/V6 , (D30) 
-a 4 (0)/\/l8 + 4(0)/\/6 (D31) 
1 4(°) a s(0) , n 6 (0) n 5 (0) 



3 Vi8 V18 
n 4 (0) 5n 3 (0) 



+ 



(D32) 



2 6 

and, where 5 assumes the values <5i and #2, 

S-(a,S) = -a\(d)/VT8 + a 6 (S)/V6, (D33) 
S+(a,d) = -a 4 {S)/Vl8 + al(6)/V6, (D34) 
1 aj(5) a 5 ((5) 5n 3 (5) 



3 # V18 6 
n 4 (6) _ n 5 (<?) n 6 (£) 



. (D35) 



Thus 



ft 



E 

6 



3G 



4(0)a 4 (<5) 



+34(<5)a 6 (0) + a\{6)a 4 (0) + 3al(0)a 6 (S)j 



+n 6 (0) - n 5 (0) - 3n 4 (0) - 5n 3 (0) 



+^(4(0K(«5) + a 5 (0)4(5) 



(D36) 



6. c at interacts with b at 6i 



Here 



S-(c) = -4(0)/Vl8 + a 6 (0)/v / 6 , (D37) 
S+(c) = -a 4 (0)/\/l8 + 4(0)/v / 6 , (D38) 
1 4(0) a B (0) n 6 (0) 



Sz(c) 



3 Vl8 718 6 
n 5 (0) n 4 (0) 5n 3 (0) 



(D39) 



and 



S_(Mi) - 24(*5i)/Vl8-2o 6 (<5i)/V6 , (D40) 
S+(Mi) = 2a 4 (5i)/\/T8-24(5i)/v / 6 , (D41) 

6 v/18 v/18 3 



(D42) 



Thus 



^=j[-4(0)a 4 (5)-34(5)a 6 (0) 
-al(S)a 4 {0) - 3oJ(0)o6(«)]/18 
+j[-4n 3 ((5) + 4n 5 (5) + 8n 6 (<*) 
-ne(0) + n B (0) + 3n 4 (0) + 5n 3 (0)]/36 
+j[-4(0)a 5 ((5) - a 5 (0)4(<5)]/9 . (D43) 

7. Summary 

Summing all the above contributions we get the 
Hamiltonian for the band at energy J for the ferro- 
conhguration as 



ft = 2fc[n 6 (R) - n 5 (R) - 3n 4 (R) - 5n 3 (R) 
r L 

-5n 3 (Ri) - 3n 4 (Ri) - n 5 (Ri) + n 6 (Ri)]/9 
+j[7n 6 (R) + 5n 5 (R) + 3n 4 (R) + n 3 (R) 
+n 3 (Ri) + 3n 4 (Ri) + 5n B (Ri) + 7n 6 (Ri)]/18 
+ ^2j[-al(R)a 4 (n + 5) -4(R + 5)a4(R) 

-34(R)ae(R + 5) - 34 (R + 5)a 6 (R) 
-24(R)a 5 (R + 5) - 24 (R + 5)a 5 (R)]/18 
+/c[4(R)a 4 (R + 5) + 4(R + <5)a 4 (R) 




FIG. 10: As Fig. 1 for a covering with one trimer per 
unit cell. In this case the paramagnetic unit cell contains 
a single trimer. The magnetic unit cell in the presence 
of antiferromagnetic trimer ordering is the same as that 
in Fig. [JJand contains two trimers. 



In Fig. QJJ we show another covering of the 
Kagome lattice with trimers. If the trimers are an- 
tiferromagnetically ordered, then the magnetic unit 
cell is the same as that of Fig. [JJ and one can ver- 
ify that the spectrum within the ground manifold is 
again given by Eqs. ([5} and (|T0|) . 
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